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Summary/Overview 

Experimental,  theoretical,  and  computational  studies  of  various  combustion  phenomena  were 
carried  out  for  a  wide  range  of  conditions.  The  parameters  that  were  considered  include  the  fuel 
type,  reactant  composition,  flame  temperature,  pressure,  and  combustion  mode.  Flames  of  FF, 
CO,  as  well  as  neat  and  mixtures  of  gaseous  and  liquid  hydrocarbons,  were  studied 
experimentally  in  the  counterflow  configuration.  The  experimental  data  were  used  as  the  basis 
for  developing  and  validating  reaction  and  molecular  transport  models. 

During  the  reporting  period,  the  following  tasks  were  completed  (numbers  correspond  to  the 
archival  publications  that  are  shown  in  Section  7.0  of  Technical  Discussion): 

1.  development  of  a  new  non-linear  extrapolation  technique  for  the  determination  of  laminar 
flame  speeds  [1]; 

2.  quantification  of  diffusion  and  kinetics  effects  on  flame  propagation  and  extinction  in 
liquid  hydrocarbon  flames  [2]; 

3.  quantification  of  the  effect  of  fuel  blending  on  the  flame  response  [3,4]; 

4.  derivation  of  an  extensive  experimental  flame  database  for  C5-C12  n-alkanc,  cyclohexane, 
and  mono-alkylated  cyclohexane  flames  [1,5]; 

5.  fundamental  kinetic  studies  of  elementary  reactions  critical  to  a  predictive  combustion 
reaction  model  [6,7]; 
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6.  refinement  of  USC  Mech  II  to  predict  Ho/CO  combustion  at  pressures  up  to  500  atm  [8]; 

7.  development  of  detailed  and  simplified  models  for  the  oxidation  of  C5-C12  //-alkanes 

[1,9,10]; 

8.  development  of  the  method  of  uncertainty  minimization  by  polynomial  chaos  expansion 
(MUM-PCE)  [11]; 

9.  gas-kinetic  theory  analysis  of  scattering  between  gas  and  molecular  clusters  with 
application  in  the  transport  properties  of  high  molecular  weight  hydrocarbons  in  dilute 
gases  [12]. 


Technical  Discussion 

1.0  Background 

The  aim  of  three-year  research  program  that  was  completed  was  to  provide  insight  into  the 
physicochemical  processes  that  control  the  burning  behavior  of  fuels  that  are  of  relevance  to 
high-speed  air-breathing  propulsion.  This  goal  was  achieved  through  combined  experimental, 
theoretical,  and  detailed  numerical  studies  of  a  number  of  fundamental  combustion  properties. 
Experimentally,  the  phenomena  of  flame  ignition,  propagation,  and  extinction  of  systematically 
chosen  fuel/oxidizer  mixtures  were  investigated.  The  studies  considered  both  gaseous  and  liquid 
hydrocarbons  in  view  of  their  importance  towards  the  accurate  description  of  the  combustion 
behavior  of  practical  fuels.  The  experimental  studies  expanded  notably  the  parameter  space  of 
existing  archival  fundamental  combustion  data.  The  detailed  modeling  of  the  experimental  data 
derived  as  part  of  this  program  and  data  available  in  the  literature,  was  performed  by  using  well- 
established  numerical  codes  and  reaction  and  transport  models  that  have  been  developed  based 
on  a  novel  approach.  More  specifically,  the  derivation  of  the  kinetics  models  was  based  on  the 
following  four  specific  objectives:  (a)  to  develop  and  test  a  suitable  array  of  mathematical  and 
computational  tools  that  can  be  used  to  quantify  the  joint  rate  parameter  uncertainty  space  and,  in 
doing  so,  to  facilitate  the  rational  design  of  reaction  models  suitable  for  any  given  target  fuel  or 
mixture  of  these  fuels;  (b)  to  examine  the  applicability  of  the  H2/CO/C14  reaction  model, 
developed  under  the  current  and  prior  AFOSR  support,  to  predict  the  phenomena  of  flame 
ignition,  propagation,  and  extinction  along  with  flame  structures,  (c)  to  examine  the  effect  of 
inelastic  collision  on  binary  diffusion  coefficients;  and  (d)  to  test  the  hypothesis  that  there  exists 
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a  critical  reaction  model  that  can  be  used  as  the  quantitative  kinetic  foundation  for  modeling  the 
combustion  of  high  molecular  weight  hydrocarbons.  The  approaches  to  meet  the  objectives 
included  quantum  chemistry  calculations,  reaction  rate  theory,  kinetic  modeling,  and  large-scale 
optimization  using  the  response  surface  method.  The  research  that  was  completed  not  only 
produced  a  comprehensive  reaction  model  suitable  for  air-breathing  propulsion  simulations  but 
also  generated  a  set  of  new  approaches  and  computational  tools  for  rational  reaction  model 
development  and  optimization.  It  was  shown  systematically  that  not  only  the  chemistry  of 
hydrocarbons  even  at  the  C2  level  is  not  well  described  but  also  that  there  are  notable 
uncertainties  associated  with  the  H2/02  chemistry,  especially  when  attempts  are  made  to  predict 
simultaneously  low-temperature  and  high-temperature  flame  phenomena.  The  kinetics  of  H2/02 
in  general  are  considered  to  be  rather  well  understood  and  are  a  very  important  subset  of  the 
kinetics  of  all  practical  hydrocarbon  fuels.  The  uncertainties  of  the  H2/02  kinetics  can  have  a 
profound  effect  on  the  simultaneous  prediction  of  ignition  and  vigorous  burning  of  H2  and 
hydrocarbons. 
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2.0  Experimental  Approach 

The  experiments  were  conducted  in  the  counterflow  configuration  as  shown  in  Figs.  1  and  2, 
which  allows  for  the  establishment  of  steady,  laminar,  planar  premixed  or  non-premixed  flames 
and  the  measurement  of  laminar  flame  speeds,  as  well  as  ignition  and  extinction  limits.  From  a 
practical  point  of  view,  these  properties  are  measures  of  the  ability  of  flames  to  be  established 
and  sustained  under  the  extreme  conditions  encountered  in  engines.  From  a  fundamental  point 
of  view,  such  data,  being  sensitive  to  the  fuel’s  oxidation  chemistry,  can  be  used  to 
compile/validate/optimize  kinetic  models. 

An  array  of  counterflow  rigs  was  used,  and  each  one  is  comprised  of: 

•  flow  system  that  is  appropriate  for  the  handling  of  gaseous  fuels,  as  well  as  liquid 
fuels  in  their  vapor  phase; 

•  several  pairs  of  variable-diameter  heated  burners  opposing  each  other; 

•  either  a  digital  particle  image  velocimetry  (DPIV)  or  a  laser  Doppler  velocimetry 
(LDV)  system  appropriate  for  velocity  measurements; 

•  several  thermocouple  systems  appropriate  for  temperature  measurements. 

One  set  of  counterflow  burners  is  housed  into  a  stainless  steel  chamber,  shown  in  Fig.  3, 
which  can  operate  from  pressures  well  below  atmospheric  up  to  10  atm  and  is  equipped  with 
state-of-the-art  remotely  controlled  functions  for  aligning  the  burners  and  achieving  ignition. 

The  vaporization  system  consists  of  a  high-precision  syringe  pump  that  injects  the  fuels  into 
a  heated  chamber  surrounded  by  a  flow  of  heated  air  or  N2.  The  flow  rates  of  all  gases  are 
measured  using  calibrated  sonic  orifices.  The  vaporization  system  has  undergone  several 
iterations  to  achieve  optimum  performance  and  reliability  for  heavy  liquid  fuels.  The 
vaporization  chamber  is  made  out  of  glass  to  prevent  fuel  cracking.  Additionally,  a  quartz 
nebulizer  with  a  flush  capillary  lapped  nozzle  was  integrated  into  the  system  in  order  to 
introduce  the  fuel  into  the  chamber  as  a  fine  aerosol.  The  details  are  shown  in  Fig.  4.  Rapid 
vaporization  is  achieved  at  relatively  low  chamber  temperatures,  allowing  steady- state 
vaporization  and  minimizing  the  chances  of  fuel  cracking.  To  prevent  fuel  condensation,  the  fuel 
delivery  line  is  wrapped  with  electrically  heated  tapes,  and  its  temperature  is  monitored  by 
properly  positioned  thermocouples.  The  burners  are  heated  with  ceramic  heating  jackets  to 
maintain  the  unburned  mixture  temperature,  Tu,  as  high  as  500  K. 
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A  thermocouple  is  used  to  monitor  the  temperature  at  the  center  of  the  burner  exit.  Before 
each  experiment,  the  entire  system  is  allowed  to  equilibrate  and  reach  steady  state,  typically 
within  two  hours. 


Non-premix ed  flame  configuration 


Single  premixed  flame  configuration 
Air  or  N2 


Fuel+ Oxidizer 


Twin  premixed  flame  configuration 


Fuel+ Oxidizer 


Fuel+ Oxidizer 


Figure  1.  Flames  that  can  be  established  in  the  opposed-jet  configuration. 


Pressure 


Figure  2.  Schematic  of  a  typical  combustion  rig. 
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Figure  3.  High-pressure  chamber  equipped  with  counterflow  burners  and  a  state  of  the  art  DPIV  system. 


Figure  4.  Schematic  and  picture  of  the  updated  liquid  fuel  vaporization  system. 

Flame  ignition  is  achieved  by  heating  air  or  N2  in  the  upper  burner  at  temperatures  well 
exceeding  1000  K.  To  heat  air  or  N2  to  temperatures  exceeding  1000  K,  the  burner  is  replaced 
by  a  quartz  tube  equipped  with  independently  controlled  helicoidal  internal  (SiC)  and  external 
(Kanthal  AI)  electric  heating  elements  that  are  capable  of  raising  the  temperature  of  air  or  N2  up 
to  1600  K  (Fotache  et  al.  1997).  A  temperature  controller/power  supply  is  used  to  operate  the 
internal  heater.  Feedback  to  the  temperature  controller  is  provided  by  a  B-type  thermocouple 
(Pt-6%  Rh/Pt-30%  Rh)  positioned  in  the  heated  air  or  N2  stream  within  the  quartz  tube  close  to 
the  surface  of  the  inner  heater.  For  a  given  reactant  configuration  and  composition,  the 
temperature  at  the  quartz  tube  exit  is  increased  and  monitored  via  another  B-type  thermocouple 
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until  flame  ignition  occurs,  at  which  point  the  tube  exit  temperature  is  defined  as  the  ignition 
temperature,  Tign.  Aspects  of  the  ignition  system  are  shown  in  Fig.  5. 


Figure  5.  Ignition  counterflow  burner  rig. 


Figure  6.  Axial  velocity  along  the  system  centerline. 

The  measurements  included  axial  flow  velocities  and  temperature.  The  axial  flow  velocities 
were  measured  along  the  stagnation  streamline  using  either  LDV  or  DPIV.  The  flow  was  seeded 
with  0.3-um  diameter  heated  silicon  oil  droplets  that  were  produced  by  a  nebulizer.  A  typical 
axial  velocity  profile  is  shown  in  Fig.  6.  The  local  strain  rate,  K,  is  defined  as  the  maximum 
absolute  value  of  the  velocity  gradient  just  upstream  of  the  flame,  as  shown  in  Fig.  6  (Law  1988). 
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The  minimum  axial  flow  velocity  just  upstream  of  the  flame  is  defined  as  the  reference  flame 
speed,  SUKf,  as  shown  also  in  Fig.  6  (Law  1988). 


3.0  Numerical  Approach 

The  experiments  were  simulated  directly  by  including  detailed  descriptions  of  chemical 
kinetics  and  molecular  transport.  Laminar  flame  speeds,  S°,  extinction  strain  rates,  Kexi,  and  Tign 
were  computed  by  utilizing  CHEMKIN-based  codes  (Kee  at  al.  1985,  1986,  1988,  1989; 
Egolfopoulos  &  Campbell  1996).  Those  codes  have  been  upgraded  by  the  Pis  to: 

•  account  for  any  type  of  combustion  mode  (Egolfopoulos  1994a;  Egolfopoulos  & 
Campbell  1996); 

•  account  for  any  type  of  boundary  condition  (Egolfopoulos  1994a;  Egolfopoulos  & 
Campbell  1996,  Egolfopoulos  et  al.  1997); 

•  account  for  thermal  radiation  (Egolfopoulos  1994b;  Zhang  &  Egolfopoulos  2000); 

•  Solve  around  singular  turning  points  at  the  states  of  ignition  and/or  extinction 
(Egolfopoulos  &  Dimotakis  1998;  Dong  et  al.  2005;  Holley  et  al.  2006,  2009); 

•  perform  mathematically  rigorous  sensitivity  analysis  of  all  dependent  parameters, 
including  S°,  Kext,  and  7jgn,  on  all  rate  constants  and  all  binary  diffusion  coefficients 
(Dong  et  al.  2005;  Holley  et  al.  2006,  2009). 

Additionally,  efficient  post-processing  codes  have  been  developed  to  perform  a  reaction  path 
analysis  that  when  combined  with  the  sensitivity  analyses,  provides  the  required  insight  into  the 
controlling  kinetic  pathways  for  each  phenomenon  that  is  investigated. 
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4.0  Results  and  Discussion 

4.1  Development  of  a  new  non-linear  extrapolation  technique  for  the  determination  of  laminar 
flame  speeds 

To  determine  S°,  SU)ref  and  K  need  to  be  measured  first,  as  shown  in  Fig.  6.  Subsequently, 
plotting  Su,ref  against  K.  S°  could  be  determined,  in  principle,  by  extrapolating  linearly  the  ,S'u.ref 
data  to  zero  stretch,  K  =  0  (e.g..  Law  1988);  however,  it  has  been  shown  (e.g.,  Tien  and  Matalon 
1991)  that  as  K  -»•  0  the  variation  of  ,SU;ref  may  not  be  linear  with  respect  to  K  because  of  the 
dilatation.  Therefore,  the  linearly  extrapolated  values  of  S°  have  been  shown  to  be  higher  than 
those  obtained  by  using  non-linear  extrapolation  (Davis  &  Law  1998)  or  by  other  more  direct 
methods  (Vagelopoulos  &  Egolfopoulos  1998). 

A  new  non-linear  extrapolation  technique  was  introduced.  Specifically,  the  variation  of  Su,ref 
with  K  was  computed  from  first  principles,  i.e.  using  a  detailed  description  of  chemical  kinetics 
and  molecular  transport  over  a  range  of  K  for  which  there  is  no  heat  loss  to  the  burner  boundary. 
In  addition,  the  S°  value  corresponding  to  K  =  0  was  computed  in  order  to  anchor  the  SUJC(-vs.-K 
correlation.  The  computed  5u>ref  and  S°  values  then  were  fitted  as  a  polynomial  function  of  K.  It 
was  shown  that  small  but  finite  uncertainties  in  rate  constants  and  reactant  transport  coefficients 
have  a  minimum  effect  on  the  shape  of  the  SUJC\-v^.-K  correlation.  Thus,  the  polynomial  function 
can  be  translated  vertically  to  fit  the  experimental  Su,ref  best  and  determine  S°.  The  sensitivity  of 
the  shape  of  the  Su,ref-A  correlation  on  kinetics  was  tested  further  by  considering  various  kinetic 
models  that  predict  reasonably  well  S°,  and  it  was  determined  that  the  shape  of  the  non-linear 
extrapolation  curves  was  affected  minimally  by  the  choice  of  the  model.  This  finding  is 
physically  reasonable.  More  specifically,  the  response  of  Su,ref  to  fluid  dynamics  is  due  to  flame 
stretch  and  dilatation.  In  opposed-jet  flames  that  are  computed  using  similar  transport 
parameters  but  two  different  kinetic  models  that  predict  similar  S°  values,  the  effects  of  stretch 
and  dilatation  is  expected  to  be  similar.  Finally,  a  computed  SLLrc(-K  correlation  from  first- 
principle  represents  the  experimental  data  better  compared  to  asymptotic  analysis  in  which  one- 
reactant,  one-step  kinetics  and  over-simplified  assumptions  for  the  transport  coefficients  are 
used. 

Similarly  to  previous  studies  invoking  non-linear  extrapolations  (e.g.,  Davis  &  Law  1998; 
Yang  et  al.  2009),  the  present  study  confirmed  that  linear  extrapolations  result  in  higher  S° 
values.  Figure  7  depicts  the  percentage  difference  between  linearly  and  non-linearly 
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extrapolated  laminar  flame  speeds  for  atmospheric  C5  -  C12  //-alkane  flames.  The  results 
illustrate  that  the  difference  is  about  5%  for  fuel-lean  mixtures  (equivalence  ratio  (p<  1.0),  but  it 
rises  to  about  20%  for  fuel-rich  mixtures  {(j)>  1.0).  Such  differences  are  rather  large  if  reaction 
kinetics  are  to  be  validated  against  S°  data.  The  increasingly  larger  discrepancies  in  fuel-rich 
mixtures  relates  to  their  sub-unity  Lewis  number  for  the  fuels  considered  herein,  given  the  large 
molecular  weight  discrepancy  between  the  fuel  and  oxygen  (e.g.,  Law  1988;  Tien  &  Matalon 
1991). 


Equivalence  Ratio^ 


Figure  7.  Average  percentage  difference  between  S°'s  of  atmospheric  C5-C12  n-alkane/air  mixtures  as  obtained 
by  linear  and  non-linear  extrapolations.  Non-linear  extrapolations  result  in  lower  S°  values. 


4.2  Molecular  transport  and  kinetics  effects  on  flame  propagation  and  extinction  in  liquid 
hydrocarbon  flames 

Lennard- Jones  12-6  potential  parameters  were  estimated  for  //-alkanes  and  /-alkenes  with 
carbon  numbers  ranging  from  5  to  14.  These  parameters,  while  somewhat  uncertain,  form  the 
basis  for  a  comprehensive  analysis  of  the  responses  of  various  /z-dodccane  flame  phenomena  to 
kinetic  and  transport  parameters.  The  logarithmic  sensitivity  coefficients  on  both  reaction  rates 
and  binary  diffusion  coefficients  were  calculated  for  laminar  flame  speeds  as  well  as  extinction 
strain  rates  for  premixed  and  nonpremixed  flames.  Results  revealed  that  the  response  of 
premixed  flames  is  generally  insensitive  to  the  fuel  diffusivity.  Coupled  with  the  lack  of 
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sensitivity  to  the  reaction  kinetics  specific  to  the  fuel  molecule,  the  computational  results  explain 
satisfactorily  the  low  sensitivity  of  the  phenomena  of  flame  propagation  and  extinction  to  the 
variation  of  the  size  of  //-alkane  fuels.  In  contrast,  the  extinction  strain  rates  of  nonpremixed 
flames  were  found  to  be  highly  sensitive  to  the  fuel  diffusivity,  in  that  larger  fuel  diffusivities 
render  the  flames  more  resistant  to  extinction.  Hence,  the  variation  of  Kext  of  nonpremixed 
flames  with  the  //-alkane  size,  as  observed  Holley  et  al.  (2007),  is  explained  satisfactorily  by  the 
variation  of  the  mass  diffusivity  of  the  fuel. 

Figure  8  illustrates  a  summary  of  the  logarithmic  sensitivity  coefficients  of  the  model 
responses  of  /z-dodecane  flames  to  key  reaction  rate  parameters  and  binary  diffusion  coefficients. 
Laminar  flame  speeds  ( S°-(p )  and  extinction  strain  rates  for  premixed  flames  (Kext-(p)  at  several 
representative  (j)  s  and  extinction  strain  rates  for  nonpremixed  flames  [Aext-fuel+nitrogen  vs. 
oxidizer]  were  considered  for  conditions  that  are  specified  in  Table  1. 

Several  observations  can  be  made  from  the  computational  results.  The  overall  influences  of 
reaction  rates  and  diffusion  fluxes  are  about  the  same  for  the  premixed  flame  responses.  The 
sensitivity  to  diffusion  originates  primarily  from  the  O2-N2  binary  diffusivity.  Kinetically,  the 
propagation  rates  of  the  premixed  flames  are  found  to  be  insensitive  to  the  chemistry  specific  to 
/z-dodecane  as  shown  here,  since  its  decomposition  in  the  preheat  zone  of  the  flame  is  too  rapid 
to  be  rate-limiting  (You  et  al.  2009).  Rather,  the  sensitivity  spectrum  includes  only  the  rate 
coefficients  found  typically  in  H2/CO  and  small  hydrocarbon  combustion,  with 
H  +  O2  — ►  O  +  OH,  CO  +  OH  — >  CO2  +  H,  and  H  +  O2  +M  — >  HO2  +M  exhibiting  the  largest 
sensitivities  as  shown  in  Fig.  8.  The  fact  that  the  computed  responses  of  premixed  flames  are 
insensitive  to  kinetic  or  diffusion  parameters  specific  to  /z-dodecane  suggests  that  the  measured 
S°u  and  KcU  of  premixed  flames  should  be  independent  of  the  nature  and  size  of  the  alkane  fuel 
molecule.  Under  similar  conditions,  the  decomposition  kinetics  of  higher  /z-alkanes  are  similar 
to  those  of  /z-dodecane.  Additionally,  for  smaller  /z-alkanes  the  larger  diffusivities  would  cause 
the  flame  propagation  rate  to  depend  even  less  on  the  fuel  diffusion  rate. 
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Figure  8.  Logarithmic  sensitivity  coefficients  on  rate  constants  and  binary  diffusion  coefficients, 
positive  values;  dark  bars;  negative  values. 


Light  bars; 


Flame 

Description 

Laminar  Flame  Sneed  (T..  =  403  K) 

S°u  1-4 

0=1.4 

s:  i.o 

0=1.0 

S°u  0.7 

0=0.7 

Extinction  strain 

rate  of  vremixed  flame 

^ex,  1-4 

Fuel/air  (0=1.4)  vs.  N2 

Kn  1  .0 

Fuel/air  (0=1.0)  vs.  N2 

^ex,  0.7 

Fuel/air  (0=0.7)  vs.  N2 

Extinction  strain 

rate  of  nonvremixed  flame 

Kext  1-30  (air) 

Fuel/  N2=l/30  vs.  air 

Kext  1-15  (air) 

Fuel/  N2=l/15  vs.  air 

Kext  1-120  (CL) 

Fuel/  N2=  1/120  vs.  02 

Kext  1-60  (02) 

Fuel/ N2=  1/60  vs.  02 

Table  1.  Conditions  of  computed  results  shown  in  Figure  8.  For  all  cases,  the  fuel  is  n-dodecane,  the  pressure  is 
atmospheric,  and  the  temperature  of  the  fuel-containing  stream  is  403  K.  The  oxidizer  is  air  for  premixed  flames, 
and  air  or  CL  for  nonpremixed  flames.  For  nonpremixed  flames  the  fuel  stream  consists  of  n-dodecane/nitrogen 
mixtures  with  varying  molar  ratios. 
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4.3  Quantification  of  the  effect  of  fuel  blending  on  the  flame  response 

5,7 s,  text’s,  and  rign’s  were  computed  laminar  flame  speeds,  as  well  as  ignition  and 
extinction  limits  of  mixtures  of  air  with  hydrogen,  carbon  monoxide,  and  C1-C4  saturated 
hydrocarbons,  were  studied  both  experimentally  and  numerically.  The  fuel  mixtures  were 
chosen  in  order  to  gain  insight  into  potential  kinetic  couplings  during  the  oxidation  of  fuel 
mixtures.  The  experiments  were  modeled  using  the  USC  Mech  II  kinetic  model  (Wang  et  al. 
2007).  It  was  determined  that  when  hydrocarbons  are  added  to  hydrogen  flames  as  additives, 
flame  ignition,  propagation,  and  extinction  are  affected  in  a  counterintuitive  manner.  More 
specifically,  by  substituting  methane  by  propane  or  //-butane  in  hydrogen  flames,  the  reactivity 
of  the  mixture  is  reduced  both  under  pre-ignition  and  vigorous  burning  conditions.  This 
behavior  stems  from  the  fact  that  propane  and  //-butane  produce  higher  amounts  of  methyl 
radicals  that  can  recombine  readily  with  atomic  hydrogen  and  reduce,  thus,  the  rate  of  the 
H  +  O2  — »  O  +  OH  branching  reaction.  The  kinetic  model  predicts  accurately  the  experimental 
data  for  flame  propagation  and  extinction  for  various  fuel  mixtures  and  pressures.  On  the  other 
hand,  it  underpredicts,  in  general,  7lgn’s. 

Additionally,  Sf  s  of  mixtures  of  air  with  80%  /z-dodccanc  +  20%  methylcyclohexane  and 
80%  /z-dodccanc  +  20%  toluene,  on  a  per  volume  basis,  were  determined  over  a  wide  range  of 
(p’s,  at  atmospheric  pressure  and  Tu  =  403  K.  The  choice  of  the  fuel  blends  was  dictated  by  their 
anticipated  compositions  in  surrogates  of  jet  fuels.  Phenomenological  analysis  revealed  that 
57 s  of  binary  fuels  mixtures  can  be  estimated  using  S°’s  and  adiabatic  flame  temperatures  of 
flames  of  the  neat  components.  The  propagation  rates  of  flames  of  the  various  binary  fuels 
blends  were  computed  using  detailed  descriptions  of  chemical  kinetics  and  molecular  transport 
and  were  found  to  be  in  good  agreement  with  the  estimations.  Although  the  initial  fuel 
consumption  pathways  and  the  resulting  intermediates  and  radicals  may  be  different  for  each 
neat  component,  the  propagation  of  flames  of  binary  fuels  was  found  to  be  mostly  sensitive  to 
the  flame  temperature  through  its  influence  on  the  main  branching  reaction  H  +  O2  — >  OH  +  O. 
Thus,  kinetic  couplings  resulting  from  the  presence  of  two  different  fuels  appear  to  have  minor 
effect  on  flame  propagation. 
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4.4  Derivation  of  an  extensive  experimental  flame  database  for  C5-C/2  n-alkane,  cyclohexane, 
and  mono-alkylated  cyclohexane  flames 

S°’  s  and  Kexl  of  premixed  C5-C12  n-alkane  flames  were  determined  over  a  wide  range  of  c/rs 
at  atmospheric  pressure  and  elevated  TVs.  S°’s  were  obtained  using  the  non-linear  extrapolation 
technique  described  in  Section  4.1.  Compared  to  linearly  extrapolated  values,  the  laminar  flame 
speeds  obtained  using  non-linear  extrapolations  were  found  to  be  1  to  4  cm/s  lower,  depending 
on  4 >•  S°’ s  of  all  77-alkane/air  mixtures  considered  in  this  investigation  are  similar  to  each  other 
and  sensitive  largely  to  the  H2/CO  and  C1-C4  hydrocarbon  kinetics.  Additionally,  the  resistance 
to  extinction  decreases  as  the  fuel  molecular  weight  increases.  Simulations  of  the  experiments 
were  performed  using  the  recently  developed  JetSurF  0.2  reaction  model  (Sirjean  et  al.  2008) 
consisting  of  194  species  and  1459  reactions.  The  experimental  S°’s  were  predicted  with  good 
accuracy  for  all  the  77-alkanes.  The  experimental  A^xt’s  were  predicted  closely  by  the  model  for 
fuel-lean  mixtures.  For  stoichiometric  and  fuel-rich  mixtures,  the  computed  Kcfls  were 
approximately  10%  lower  than  the  experimental  values.  Insights  into  the  physical  and  chemical 
processes  that  control  the  response  of  n-alkane  flames  were  provided  through  detailed  sensitivity 
analyses  on  both  reaction  rates  and  binary  diffusion  coefficients.  Figures  9  and  10  depict  the 
experimental  and  computed  S°’s. 

S°  ’  s  were  measured  in  the  counterflow  configuration  for  cyclohexane/air, 
methylcyclohexane/air,  ethylcyclohexane/air,  77-propylcyclohexane/air,  and  n- 
butylcyclohexane/air  mixtures  at  atmospheric  pressure,  Tu  =  353  K,  and  for  a  wide  range  of  r/rs. 
Cyclohexane/air  flames  propagate  somewhat  faster  than  mono- alkylated  cyclohexane/air  flames. 
Flames  of  mono-alkylated  cyclohexane  compounds,  i.e.  methyl-,  ethyl-,  77-propyl-,  and  77-butyl- 
cyclohexane,  were  found  to  have  similar  S°’ s,  suggesting  that  the  different  alkyl  groups  have  a 
secondary  effect  on  flame  propagation.  The  experimental  data  were  modeled  using  JetSurF  1.1 
(Sirjean  et  al.  2009),  a  detailed  reaction  model  for  the  combustion  of  cyclohexane  and  its 
derivatives,  and  satisfactory  agreements  were  found.  Based  on  the  analysis  of  the  modeling 
results,  the  somewhat  lower  rates  of  mono-alkylated  cyclohexane  flame  propagation  were 
attributed  to  the  greater  production  of  propene  and  allyl  and  the  increased  H-atom  scavenging  by 
these  C3  intermediates.  Although  these  fuel-specific  reaction  kinetic  features  do  not  limit  the 
overall  oxidation  rates,  the  distribution  of  the  cracked  products  exert  influences  on  flame 


Laminar  Flame  Speed,  cm/s  Laminar  Flame  Speed,  cm/s 
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propagation,  leading  to  the  subtle  differences  in  S°’s  that  were  observed  experimentally.  Figure 
11  depicts  the  experimental  and  computed  S°’ s. 


Figure  9.  Experimental  and  computed  S°'s  at  Tu  =  353  K.  (a)  u-CLH  12,  (b)  «-C6Hi4,  (c)  n-C7H 16,  and  (d)  n- 
CLH  i 8/air  mixtures.  Symbols:  present  experimental  data  obtained  with  non-linear  extrapolation;  lines: 
simulation  results  using  the  JetSurF  0.2  kinetic  model. 
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(a)  (b) 

Figure  10.  Experimental  and  computed  S°’s  at  Tu  =  403  K.  (a)  n-CgEEo/air  mixtures  (•);  (b)  n-CioEWair 
mixtures,  (•)  present  experimental  data  obtained  with  non-linear  extrapolation;  (□)  Experimental  data  of 
Kumar  &  Sung  (2007);  lines:  simulation  results  using  the  JetSurF  0.2  kinetic  model.  The  error  bars  indicate 
2- a  standard  deviations. 


Equivalence  Ratio,  <p 


Figure  11.  Experimental  and  computed  S°’s  of  methyl-,  ethyl,  n -propyl,  and  n-butyl-cyclohexane/air  flames  at 
p  -  I  atm.  Symbols:  (■)  present  experimental  data  at  Tu  =  353  K,  (A)  Dubois  et  al.  (2009).  Lines:  ( — )  simulation 
using  Model  I  (JetSurF  1.1)  at  Tu  =  353  K,  (— )  simulation  using  Model  I  (JetSurF  1.1)  at  Tu  =  403  K.  The  error  bars 
indicate  2- a  standard  deviations. 
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4.5  Fundamental  kinetic  studies  of  elementary  reactions  critical  to  a  predictive  combustion 
reaction  model 

The  kinetics  of  the  reaction 

CO  +  H02*  -►  C02  +  *OH  (Rl) 

was  studied  using  a  combination  of  ab  initio  electronic  structure  theory,  transition  state  theory, 
and  master  equation  modelling.  The  potential  energy  surface  was  examined  with  the  coupled 
cluster  CCSD(T)  method  and  complete  active  space  CASPT2  methods  (Andersson  et  al.  1990, 
1992).  The  classical  energy  barriers  were  computed,  as  shown  in  Fig.  12.  The  energy  barriers 
were  found  to  be  18  and  19  kcal/mol  for  CO  +  HCF*  addition  following  the  irons  and  cis  paths, 
respectively. 


rS6  35.1  FS7  35.7 


Figure  12.  Potential  energy  diagram  for  CO  +  H02*  — >  products.  For  CO  +  H02*  — >  CO 2  +  *OH,  the  energy 
values  were  determined  using  the  CCSD(T)/CBS  method,  and  included  zero-point  energy  corrections.  For 
CO  +  H02*  — *  HC*0  +  02,  the  energy  values  were  taken  from  Martlnez-Avila  et  al  (2003). 

For  the  cis  path,  rate  constant  calculations  were  carried  out  with  canonical  transition  state  theory. 
For  the  trans  path,  master  equation  modelling  also  was  employed  to  examine  the  pressure 
dependence.  Special  attention  was  paid  to  the  hindered  internal  rotations  of  the  HOOOO  adduct 
and  transition  states.  The  theoretical  analysis  showed  that  the  overall  rate  coefficient  is 
independent  of  pressure  up  to  500  atm.  The  computed  rate  coefficient  was  found  to  be 
substantially  smaller  than  some  of  the  rate  recommendations  in  the  literature,  as  shown  in  Fig. 
13.  Based  on  the  current  theoretical  calculations,  the  following  rate  expression  was 
recommended  for  reaction  Rl 

kx ( cm 3 /mol  -  s)  =  1 ,57xl05  r21  V9030^ 
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for  300  <  T  <  2500  K  with  the  uncertainty  factor  equal  to  8,  2  and  1.7  at  temperatures  of  300, 
1000,  and  2000  K,  respectively. 


Figure  13.  Rate  coefficient  of  reaction  Rl. 

Q 

The  gas-phase  reaction  of  benzene  with  OOP)  is  of  considerable  interest  for  modeling  of 
aromatic  oxidation  and  because  there  exist  fundamental  questions  concerning  the  prominence  of 
intersystem  crossing  in  the  reaction.  While  overall  rate  constant  of  reaction  (Rl)  has  been 
studied  extensively,  still  there  are  significant  uncertainties  in  the  product  distribution.  The 
reaction  proceeds  mainly  through  the  addition  of  the  O  atom  to  benzene,  forming  an  initial  triplet 
diradical  adduct  that  either  can  dissociate  to  form  the  phenoxy  radical  and  H  atom  or  undergo 
intersystem  crossing  onto  a  singlet  surface  followed  by  a  multiplicity  of  internal  isomerizations 
leading  to  several  possible  reaction  products.  In  this  work  the  product  branching  ratios  of  the 
reaction  were  studied  over  the  temperature  range  of  300  to  1000  K  and  pressure  range  of  1  to 
10  Torr.  The  reactions  were  initiated  by  pulsed  laser  photolysis  of  NO2  in  the  presence  of 
benzene  and  helium  buffers  in  a  slow-flow  reactor,  and  reaction  products  were  identified  by 
using  the  multiplexed  chemical  kinetics  photoionization  mass  spectrometer  operating  at  the 
Advanced  Light  Source  (ALS)  of  the  Lawrence  Berkeley  National  Laboratory.  Phenol  and 
phenoxy  radicals  were  detected  and  quantified.  Cyclopentadiene  and  cyclopentadienyl  radical 
were  identified  directly  for  the  first  time.  Finally,  ab  initio  calculations  and  master 
equation/RRKM  modeling  (Holbrook  et  al.,  1996;  Gilbert  and  Smith  1990)  were  used  to 
reproduce  the  experimental  branching  ratios,  yielding  pressure-dependent  rate  expressions  for 
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the  reaction  channels,  including  phenoxy  +  H,  phenol,  and  cyclopentadiene  +  CO,  as  shown  in 
Fig.  14.  These  computed  rate  coefficients  were  used  to  update  the  USC_Mech  II  reaction  model 
(Wang  et  al.,  2007)  so  that  predictions  of  aromatics  combustion  are  improved. 
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Figure  14.  Arrhenius  plots  for  the  reaction  channels  leading  to  (a)  phenol,  (b)  phenoxy  +  H,  (c)  cyclopentadiene  + 
CO,  (d)  benzene  oxide,  (e)  2,4-cyclohexadienone,  and  (f)  2,5-cyclohexadienone.  Symbols  are  rates  computed  with 
Monte  Carlo  RRKM/master  equation  modeling  and  lines  represent  fitted  Arrhenius  expressions. 
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4.6  Refinement  of  USC  Mech  II  to  predict  H2/CO  combustion  at  pressures  up  to  500  atm 

The  high-pressure  oxidation  of  inert-diluted  CO/O2  mixtures  doped  with  150-200  ppm  of  FT 
was  studied  behind  reflected  shock  waves  in  the  University  of  Illinois  at  Chicago  high-pressure 
single  pulse  shock  tube.  The  experiments  were  performed  for  stoichiometric  (0=1)  and  fuel 
lean  (0  =  0.5)  mixtures  over  the  temperature  and  pressure  ranges  of  1000-1500  K  and  21- 
500  bar,  respectively.  Stable  species  sampled  from  the  shock  tube  were  analyzed  by  gas 
chromatographic  (GC)  and  GC/mass  spectrometric  techniques.  The  experimental  data  were 
simulated  using  the  USC  Mech  II  kinetic  model  (Wang  et  al.  2007).  These  simulations  showed 
that,  within  the  experimental  uncertainty,  USC  Mech  II  was  able  to  capture  the  experimental 
trends  for  the  lower  pressure  data  sets  (average  nominal  pressures  of  24  and  43  bars);  however 
the  model  underpredicted  the  CO  and  O2  decay  and  subsequent  CO2  formation  for  the  higher- 
pressure  data  sets  (average  nominal  pressures  of  256  and  450  bars).  The  elevated  pressure  data 
sets  span  a  previously  unmapped  regime  and  have  served  to  probe  HO2  radical  reactions,  which 
appear  to  be  among  the  most  critical  reactions  under  the  experimental  conditions  that  were 
considered.  With  updated  rate  parameters  for  a  key  HO2  radical  reaction,  USC  Mech  II  was 
shown  to  reconcile  the  elevated  pressure  data  sets,  thereby  extending  the  model  predictability  to 
an  extreme  range  of  conditions. 


400 

300 

200 


C 

.2 


£  500 


a  400 
S  300 
200 
100 
0 


24  atm.  <S»=1  43  atm.  <4=1 


d,  A  .«ri 

“  T  *4  11  mrr  >. 

B  a 

400 

300 

200 

100 

0 

5* 

PH  rcfl 

256  atm.  4>=I 

450  atm.  4=1 

1®  ij* 

400 

300 

*  ♦  St  *  Jeair, 

H  B  B  B 

200 

100 

0 

"s® D 

B  B  H  H  S*5 

1100  1200  1300  1  400  1500  1100  1200  1300  1400 

Temperature/  K 


Figure  15.  Experimental  data  for  CO  (•)  and  CO2  (■)  measured  by  shock  heating  of  a  mixture  containing  500  ppm 
CO,  200  ppm  H2,  02  and  argon  with  the  equivalence  ratio  equal  to  1.0.  The  computational  results  are  □:  base  USC 
Mech  II,  x;  updated  USC  Mech  II. 
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4.7  Development  of  detailed  and  simplified  models  for  the  oxidation  of  C5-C12  n-alkanes 

A  detailed  kinetic  model  was  proposed  for  the  combustion  of  n-alkanes  up  to  /z-dodecane  and 
for  temperatures  above  850  K.  The  model  was  validated  against  experimental  data,  including 
fuel  pyrolysis  in  plug  flow  and  jet- stirred  reactors,  laminar  flame  speeds,  and  ignition  delay  times 
behind  reflected  shock  waves,  with  n-dodecane  being  the  emphasis.  Selected  model  validation 
results  are  shown  in  Figs.  16  and  17.  The  model  was  used  as  the  basis  for  JetSurF  0.1  and  1.0 
with  a  wider  range  of  validation  tests  reported  on  the  JetSurF  web  site  (Sirjean  et  al.  2008,  2009) 
and  in  Figs.  9  and  10. 


Residence  Time  (s) 


Figure  16.  Experimental  (symbols:  Herbinet  et  al.  2007) 
and  computed  (solid  line:  detailed  model;  dashed  line: 
simplified  model)  conversions  of  /7-C12H26  as  a  function 
of  the  residence  time  for  the  pyrolysis  of  2%  77-C12H26  in 
He  in  a  jet-stirred  reactor  (p  =  1  atm). 


1000  K  /  7 

Figure  17.  Experimental  (symbols)  and  computed  (solid 
line:  detailed  model;  dashed  line:  simplified  model  not 
shown  here)  ignition  delay  times  for  H-C10H22  oxidation 
behind  reflected  shock  waves. 


Analysis  of  the  computational  results  revealed  that  for  a  wide  range  of  combustion 
conditions,  the  kinetics  of  fuel  cracking  to  form  smaller  molecular  fragments  is  fast  and  may  be 
decoupled  from  the  oxidation  kinetics  of  the  fragments.  Subsequently,  a  simplified  model 
containing  a  minimal  set  of  4  species  and  20  reaction  steps  was  developed  to  predict  the  fuel 
pyrolysis  rate  and  product  distribution.  Combined  with  the  base  C1-C4  model,  the  simplified 
model  predicts  fuel  pyrolysis  rate  and  product  distribution,  laminar  flame  speeds,  and  ignition 
delays  as  accurately  as  the  detailed  reaction  model. 
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4.8  Development  of  the  method  of  uncertainty  minimization  by  polynomial  chaos  expansion 

The  ability  of  a  reaction  model  to  predict  the  combustion  behavior  of  a  fuel  relies  on  the 
rigorous  quantification  of  the  kinetic  rate  parameter  uncertainty.  Although  the  accuracy  of  a 
detailed  kinetic  model  may  be  ensured,  in  principle,  by  a  multi-parameter  optimization,  the 
inherent  uncertainties  in  the  fundamental  combustion  targets  used  for  optimization  cause  the 
derived  optimized  model  to  be  characterized  by  a  finite  kinetic  parameter  space.  In  this  work, 
spectral  expansion  techniques  were  developed  and  employed  to  quantify  these  uncertainties.  The 
resulting  method,  called  the  Method  of  Uncertainty  Minimization  by  Polynomial  Chaos 
Expansion  (MUM-PCE)  was  used  to  analyze  a  detailed,  H2/CO/C1-C4  kinetic  model  for  ethylene 
combustion  (Wang  et  al.,  2007).  Uncertainty  was  quantified  for  both  the  as-compiled  model  and 
the  optimized  model  and  was  propagated  subsequently  into  a  wide  variety  of  combustion 
experimental  conditions.  Examples  for  laminar  flame  speeds  and  shock  tube  ignition  delays  are 
shown  in  Figs.  18  and  19. 


Figure  18.  Experimental  data  (o;  Egolfopoulos  et  al. 
(1990),  •:  Jomaas  et  al.  (2005),  0:  Hassan  et  al.  (1998)) 
and  2a  uncertainty  bands  computed  for  the  laminar 
flame  speed  of  ethylene-air  mixtures  at  p  =  1  atm. 


1 0,000/7  (1/K) 

Figure  19.  Experimental  data  (symbols.  Brown  and 
Thomas,  1999)  and  computed  (dashed  line:  as-compiled, 
unoptimized  model;  solid  line:  optimized  model,  bands;  2a 
model  uncertainties)  ignition  delay  time  of  l%C2H4-3%02- 
Ar  mixtures  behind  reflected  shock  waves. 
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4.9  Gas-kinetic  theory  analysis  of  scattering  between  gas  and  small  particles  with  application  in 
the  transport  properties  of  high  molecular  weight  hydrocarbons  in  dilute  gases 

A  theoretical  study  was  conducted  to  examine  the  uncertainties  in  the  transport  theory  of 
molecular  clusters  and  small  particles  in  laminar  flow  regime  (Li  and  Wang  2003a,  2003b, 
2005).  The  theory  features  a  rigorous  gas  kinetic  theory  analysis,  considers  the  effect  of  non- 
rigid  body  collision,  and  has  been  shown  to  reproduce  the  Chapman-Enskog  theory  (Chapman 
and  Cowling,  1970)  of  molecular  transport  in  the  small  particle  size  limit,  Epstein's  model 
(Epstein,  1924)  of  particle  drag  in  the  rigid-body  limit,  and  the  Stokes-Cunningham  equation 
(Cunningham,  1910)  for  the  drag  of  micrometer  size  particles.  A  key  uncertainty  in  the 
application  of  the  generalized  theory  lies  in  the  uncertainty  of  momentum  accommodation  of 
molecules  or  particles  in  the  size  range  of  a  few  nanometers.  The  cause  for  this  uncertainty  was 
attributed  to  a  transition  in  the  dominant  mode  of  gas  particle  interactions  from  specular-like 
scattering  to  diffuse  scattering.  As  an  example,  Fig.  20  depicts  this  transition  by  comparing  the 
theoretical  reduced  collision  integrals  for  polyethylene  glycol  (PEG)-air  molecule  interactions 
and  the  experimental  ones  derived  by  Fernandez  de  la  Mora  and  coworkers  (Nasibulin  et  al., 
2002;  Ude  and  Fernandez  de  la  Mora,  2003).  In  order  to  apply  accurately  a  generalized  transport 
theory  for  large  molecules,  a  substantially  larger  set  of  experimental  mobility  and  diffusivity  has 
to  be  available. 


R  (nm) 

Figure  20.  Variations  of  the  reduced  collision  integrals  as  functions  of  the  radius  of  PEG  nanoparticles.  Symbols 
are  experimental  data  observed  in  air  at  room  temperature  (diamonds:  Nasibulin  et  al.,  2002;  circles:  Ude  and 
Fernandez  de  la  Mora,  2003). 
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Technical  Interactions  with  AFRL  Researchers 

Several  meetings  with  Dr.  Tim  Edwards  of  AFRL  took  place  during  the  reporting  period. 
The  purpose  of  these  meetings  was  to  discuss  the  measurements  on  the  propagation,  ignition,  and 
extinction  limits  of  flames  of  neat  hydrocarbon  fuels  as  well  as  a  variety  of  petroleum-derived 
and  synthetic  jet  fuels  that  were  supplied  to  the  USC  Combustion  and  Fuels  Laboratory  by  Dr. 
Edwards.  The  discussions  has  focused  on  both  the  experimental  and  modeling  work  pertaining 
to  jet  fuels  and  their  surrogates,  as  well  as  ongoing  work  on  smaller  hydrocarbons  in  the  C1-C4 
range. 


Other  Technical  Interactions 

A  close  interaction  was  established  with  Professors  Ron  Hanson,  Tom  Bowman  and  Heinz 
Pitsch  of  Stanford  University,  Professor  C.K.  Law  of  Princeton  University,  Professors  Nick 
Cemansky  and  David  Miller  of  Drexel  University,  Professor  Angela  Violi  of  University  of 
Michigan,  Dr.  Wing  Tsang  of  National  Institute  of  Standards  and  Technology,  and  Dr.  Meredith 
Colket  of  United  Technologies  Research  Center  on  the  experimental  and  modeling  aspects  of 
this  research.  Additionally,  collaborative  research  efforts  were  initiated  with  Dr.  Stephen 
Klippenstein  of  the  Argonne  National  Laboratory  and  Dr.  Craig  Taatjes  of  the  Sandia  National 
Laboratories 
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Inventions,  Technology  Transitions  and  Transfers 

Part  of  the  work  described  in  this  report  resulted  in  the  following  technology  transition  and 
transfer: 


Performer:  Dr.  Hai  Wang,  University  of  Southern  California,  (213)  740-0499 


Customer: 

Phone: 

Result: 

Application: 


Dr.  Meredith  Colket,  United  Technology  Research  Center,  411  Silver  Lane 
East  Hartford,  CT  06108 
(860)  610-7000 

Reaction  model  for  predicting  aromatics  and  soot  formation  in  fuel  rich  flames 
Low-emission  gas  turbine  engines 


